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Abstract 

We calculate the hadronic coupling constants §nn-k an d g p „ in QCD, including dynamical 
quarks in the framework of staggered fermions in the lattice approach. For the nucleon- 
pion coupling we obtain §nn-k = 13.8 ± 5.8, to be compared with the experimental value 
13.13 ± 0.07 [1]. The pirir coupling has been analysed for two different sets of operators 
with the averaged result g p „ = 4.2 ±1.9 which is to be compared with the experimental 
value 6.06 ± 0.01 [2]. 
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1 Hadronic coupling constants in lattice QCD 



Recent lattice QCD calculations have given encouraging results for the hadronic mass 
spectrum. This success motivates the launching of more complicated analyses, in particular 
the calculation of dynamical properties of hadrons. The aim of the present analysis is 
the calculation of coupling constants which parameterise the strength of the interaction 
between hadrons at low energies, above all the classical nucleon-pion coupling constant. 

We have based the analysis on the staggered fermion action with four degenerate 
flavours of dynamical fermions. The lattice size is L 3 X T = 16 3 X 24 and we have generated 
85 configurations at (3 = 5.35, m = 0.01 using the hybrid Monte Carlo algorithm. Details 
with respect to the production of configurations can be recalled from ref.[3]. 

Hadronic coupling constants can be determined by calculating 2- and 3-point functions 
on the lattice using the path integral method. Compared to the analysis of the mass 
spectrum, however, the calculation of coupling constants is, from the theoretical as well 
as from the numerical point of view, much more involved. 

Above all we have to choose a set of operators very carefully so that the following 
criteria are fullfilled: First the operators should have a good overlap with the hadronic 
particle state under consideration and, if possible, a high suppression of the unwanted par- 
ity partner state created by the usual quasi-local staggered operators. Secondly, because 
the continuum limit of the staggered lattice action corresponds to a theory describing four 
mass degenerate quark flavours the operators have to be chosen such that the coupling is 
allowed by SU{A) F flavour symmetry group selection rules. 

Next one has to find a function to parameterise the 3-point function in which the 
hadronic coupling constant appears as a parameter. Besides the usual transfer matrix 
techniques and the LSZ formalism, which relates the correlation functions to coupling 
constants, one has to determine the corresponding SU{A) F Clebsch-Gordan (CG) coeffi- 
cients because the coupling constants are defined as flavour independent numbers. This, 
however, requires in a first step to relate the lattice operators to the corresponding oper- 
ators in the continuum so that SU{A) F quantum numbers can be assigned to the states. 

Due to the technically complicated analysis, we proceed in two steps and consider the 
g pjm and gNNw coupling constants separately (see also ref.[4] and [5]). 

2 The g pmr coupling constant 

The effective interaction Lagrangian which parameterises the decay of a p meson into two 
pions is of the form 



Here c abc is the SU(4:)f Clebsch-Gordan coefficient which ensures that the Lagrangian be- 
haves like a scalar under flavour transformations; g p „ does not depend on flavour quantum 
numbers. 

For the p meson we used the local operator [No.4] 1 because it delivered quite stable 
results in the analysis of the mass spectrum . Moreover, the numerical evaluation of the 3- 

J The numbers of the operators refer to the numbers used in ref.[3]. 
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point function is much eased by using local operators. However, like most of the staggered 
hadron operators the p operator has the unwanted property of coupling also to a state 
with opposite parity, namely a b\ meson. This parity-partner state contributes to the 
correlation functions so that additional terms have to be taken into account, thus leading 
to more complicated parameterisations. Therefore we used also a non-local-in-time (NLT) 
p operator [No. 23] which belongs to the same representation of the lattice symmetry group 
as the local p but whose coupling to the parity partner meson is dynamically suppressed. 

One of the pion operators (ir b ) was chosen to be the lattice Goldstone (LG) pion [No. 2] 
which, besides locality, has the additional advantage of being the lightest pionic state in 
our simulation 2 so that the kinematics of the p decay is closest to the situation in the 
continuum. Moreover, the parity partner of the pion created by this operator corresponds 
to an exotic state and is highly suppressed so that it can be neglected in the analysis. 

With these operators given, the third operator (ir c ) has to be chosen in a way that 
the CG coefficient c abc is non-zero. Therefore we adopted a non-local pion, the 1-link 
local-in-time (LT) operator [No. 7] which has an exotic, suppressed parity partner. We 
then calculate the 3-point correlation function 

r 3 (*i, * 2 ) = ^{qpWatfp, , ^K(O)), (2) 

being the polarisation vector of the p meson. The tilde indicates that the corresponding 
field operator is defined in momentum space, e.g. 

^-'l^E^V), (3) 

xev 

where V is the spatial volume of the lattice. 

The time slice t 2 of the pion has been fixed to t 2 = 4 and we evaluated the correlation 
function for all times t x of the p meson. The p meson was given the minimal lattice 
momentum q p = ^f-e z in the z direction, and q w was set to 0. The pion source 7r then 
projects on a state with non- vanishing momentum. 

To define a parameterisation of this correlation function we insert complete sets of 
states in eq.(2). Assuming large time distances between the operators we can restrict 
ourselves to the low lying particle states and the following expression can then be derived 
for T > h > t 2 > 0: 

T 3 (t u t 2 ) « (^,^1^(0)1^, ?;)y^Ae- £ - (il -' 2) e-^ i2 (-l)' 1+ H 

(7r c ,^|^(g,)|7r i ,0)y^X e -^( tl - t =) e - £ -( T - tl )(-l) tl+t2 + 
(7r i ,0|7r c (g,)| / o a ,g,)^Z, a e- B -( T - tl )e-^ ta (-l) tl+ta + 

(n, o|f c (g;)|6 la ,g;)^z ilae -^«( T - ti ) e - £ - t2 (-i) t2 + 

<&i a ,?„|7r 6 (0)|7r c , q p )^Z bla Z Wc e- E ^^-^e- E '^. (4) 

2 This is due to the fact that our lattice is not close enough to the continuum to completely restore the 
flavour symmetry in the chiral sector. 
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Z M denotes the 1-meson amplitude of the operator Ai, \J Z M = {0\A4\ Meson ). These 
constants together with the masses and energies of the different particle states are taken 
from the analysis of 2-point functions. The corresponding formula for T t 2 t\ 
is 

T 3 (t u t 2 ) « (7r i ,0| / 5 a (g,)|7r c ,g,)^ c Z Wl e-^( ta - tl )e- B - tl (-l) tl+ta + 
(p a , q P \n c (q p )\TT b , 0}^/Z Wh Z Pa e- E ^e- E ^ T -^(-iy^ + 
(7r c ,g^ i (0)| / > a ,g;)yZ^ e - £ -^- tl ) e - £ -( T - t2 )(-l) tl+t2 + 
(b la ,q p \t c (q p )\ir b , 0)^Z ila e"^- tl e-^( T - ta )(-l) ta + 
(7r c , ^1^(0)16^, g^^.Z^e-^-^^e-^^). (5) 

In principle eq.(4) should include terms accounting for the parity partners of the 
mesons. However, as the opposite parity states to both pions are exotic and therefore 
suppressed (i.e. Z ~ 0, as can be seen in the analysis of 2-point functions), they need not 
be taken into account. On the other hand, the coupling described by the matrix element 
where the p meson parity partner and the two pion states are involved, violates lattice 
symmetry and is therefore forbidden. For the NLT p also the last two terms in eq.(4) can 
be neglected as the the coupling of this operator to the b\ state is suppressed. 

In order to find the relation between the matrix elements, e.g. (p a |7Tj|7r c ), and the cou- 
pling constants we apply the LSZ reduction formula and evaluate the resulting Green's 
function according to the effective interaction Lagrangian in eq.(l). The resulting expres- 
sion consists of a free propagator for the ir b particle, a kinematic factor K(q p ) for the 
truncated vertex and the hadronic coupling constant: 



(p a ,q P \K b (0)\ir c ,q P ) = g pm c a bcK(q P )\/ Z nb P((E Pa - E n J,E nb ). (6) 

P(E,uj) denotes the free lattice particle propagator for a particle with pole energy u 
propagating at energy E which is defined through 

P(E,u)= / nh " ■ (7) 
coshw — cosh E 

The kinematic factor for the process p kit \s given by 



R \1p) = \1p\— 

" l Pa 

The parameter function was fitted to the numerical data by a procedure minimizing 
X 2 - The masses and energies of the various states could be kept fixed as they are already 
known from the analysis of 2-point functions. Thus we have one fit parameter for the 
correlation function with the NLT p operator and 3 fit parameters when using the LT p 
operator. However, in the latter case only one parameter, namely g p „ is of physical interest 
whereas the others are proportional to the matrix elements where the pion operator serves 
as an interpolating field between a b\ and a 7r state. The time interval was varied in the 
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fit procedure; however the variation could be done only in a restricted range since due to 
the involved parameter function the fit soon becomes unstable when the number of fitting 
points is reduced. The result of the fit is a product of the coupling constant g p7m , the 
amplitudes Z and the CG coefficients. The Z factors, like the masses, are known from 
2-point functions and can readily be divided out. 

To derive the value for g p7m , we still have to determine the CG coefficients c abc . The 
flavour non-singlet meson states are members of an SU(4:)f 15-plet so that we need to 
calculate CG coefficients for the product of representations 15 ® 15 — ► 15. The CG 
coefficients have already been calculated for this case in ref.[6], but we still have to relate 
the lattice operators to continuum operators. To do so, we write the staggered mesonic 
operators in the following form 

M SF (y) = q{y)[T s ® T F ]q(y), (9) 

where T s , T F define the spin and flavour content of the meson and q denotes a staggered 
quark field defined on a hypercube of the lattice: 

q aa (y) = J2, T v a x^y + v)- (io) 

T^ a = (ji 1 ■ ■ -74 4 ) aa is a product of 4x47 matrices and a, a denote spin and flavour indices 
of the quark, respectively. To assign continuum flavour values to the lattice operators we 
call the a = 1 component of the lattice quark operator a u quark, the a = 2 component 
a d quark, the a = 3 component an s quark and the a = 4 a c quark. 3 If we express the 
flavour matrix T F in eq.(9) through SU(4:)f Gell-Mann matrices, we are able to identify 
the lattice meson operators with continuum operators used in ref.[6]. On the other hand, 
the relationship between the lattice meson operator in eq.(9) and the lattice operators 
given in ref.[3] can be readily found by using the above formulae and trace theorems for 7 
matrices. 

As the T F matrices correspond to linear combinations of SU(4:)f Gell-Mann matrices 
the staggered meson operators correspond to linear combinations of continuum operators 

K att ~ £(.fM)X Cont - (11) 

b 

M^ att denotes a member of the lattice meson multiplet and M b Cont is a member of the 
continuum multiplet. Thus the lattice coefficient c abc in eq.(6) is given by 

c aic = Y,(fMfMf«ycCaeT, (12) 

d.ej 

where the continuum CG coefficients can be taken from ref.[6]. To be specific, we give 
the correspondence of the lattice meson operators, relevant in the evaluation of the g p „ 

3 This assignment is arbitrary as all staggered quark flavours are mass degenerate. However, consistent 
definitions have to be used when group theoretical relations including different operators are evaluated. 
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coupling constant, to the continuum meson operators: 

pLatt ^ l(K* + - K*~ + D* + - D*~) (13) 

irf att <-> ^(K+ + K~ - D+ - D~). 

These relations determine the matrices (/m)„ an( i thus allow the derivation of the CG 
coefficient c abc , 

Cabc = V2. (14) 

Now we are ready to determine the g p „ coupling constant from the numerical data 
which are shown in fig.l for the case of the local p operator. We find for the two different 
cases of the local and the NLT p operator 

g pim = 3.6 ±2.5 (LT), 

g pm = 4.7 ±2.8 (NLT), (15) 

from which we derive the average value 

^„=4.2±1.9. (16) 

Compared to the experimental value of g p „ = 6.06 ± 0.01 our results are too small but 
still compatible with the measurement within la of the large errors. Moreover, our results 
are larger than the results in ref.[7] based on Wilson fermions for a smaller lattice which 
spread between g p „ ~ 2.8 and g p „ ~ 3. 

3 The qnn-k coupling constant 

The effective Lagrangian associated with the NNir coupling is defined as 

= gNN.c' abc N a (x) l5 N b (x)^(x). (17) 

This equation holds only for SU(2) F flavour symmetry subgroup or, more precisely, if 
the nucleon belongs to a 2-dimensional representation and the meson to a 3-dimensional 
representation of the same SU(2) F subgroup of the full flavour symmetry group. For 
the full SU(4:)f flavour symmetry in a general case, however, we have to consider two 
independent coupling constants gi and g 2 so that eq.(17) must be interpreted as 

m f t (x) = ( gi cl bc +g 2 cl bc )N a (x) l5 N\x)7r%x). (18) 

Mathematically, this is due to the fact that the product of the 20-dimensional representa- 
tion to which the baryon belongs in SU(4:)f , and the 15-dimensional meson representation 
reduces, among others, to two different 20-dimensional baryon representations. This situ- 
ation is analogous to SU(3)f where F and D couplings, gp~ and a g£, must properly be 
distinguished. 
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The set of operators we chose in this case consists of a local polarised baryon opera- 
tor, the lattice Goldstone pion [No. 2] and a baryonic wall source operator. The baryon 
should be polarised since the contribution of the NNir coupling to the correlation function 
vanishes otherwise. The Goldstone pion again was used because of its low mass, its local 
character and the suppressed opposite parity state. 

We decided to use a wall source operator for the second baryon as it has a much 
better overlap with the ground state particle compared to the local operator. This is 
especially important for the case of 3-point functions because time distances between the 
operators are of course even smaller if three operators are defined on the lattice so that 
the suppression of higher states by the exponential functions is worse. Indeed, the quality 
of the numerical data improved significantly compared to earlier attempts where we used 
a local source. Moreover, the results of the analysis of 2-point functions, especially the 
amplitude Z N , are much more stable when the wall source operator is employed. 

The 3-point function we evaluated numerically, is of the form 



where N denotes a nucleon operator with positive polarisation in the z direction and N w 
is the nucleon wall source operator. Since the wall source operator we used projects onto 
states with vanishing momentum we have to give opposite momenta to the nucleon sink 
and the meson operator (the NNir coupling vanishes as p — ► 0). p was chosen to be the 
minimal lattice momentum in the z direction. The time of the sink nucleon was fixed at 
t 1 = 8 and the correlation function has been evaluated for all time slices t 2 of the pion. 

A parameterisation of this function can again be found by using transfer matrix tech- 
niques and inserting a complete set of particle states between the operators. Neglecting, 
for large time distances between the operators, the influence of higher mass states, we find 
the following parameterisation for T t 2 t\ ^0: 




(19) 





( ir b ,-p\ N a (p) | JV c ,0 ) ^Z %h Z Nc e- E ^-^e- E ^(-lf^ + 



( N a ,p\ NY(S) | ir b ,p) ^Z Wb Z Na e- E ^ T -^e- E ^(-iyTc + 



( N a ,p\w b (-p) \N e ,6) ^Z Na Z Nc e- E ^ T -^h- E ^-^\-l)^Cc^c + 



( n ,-p\ N a (p) \N>J) ^fZ Wb Z N ,e- E ^-^e- E ^(-iy* + 



( N' a ,p\lf b (-p) | JV' c ,0 ) ^Z K Z N ,e- E K (T - t2) e- E < (t2 - tl '>Cc^c + 



( N' e ,6\ M-p) | N a ,p) yvz^ e -^ (T - t2) e -^( t - ti )(-i)( t2+ti )es'cT C + 



( JV c ,0 | W b (-p) | N' a ,p) ^Z NL Z Nc e- E ^ T -^e- E K (t ^\-l)^Cc<rc, 
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and for T > t 1 > t 2 > 0: 



(21) 



N a ,p x b (-p) 



N„0 



Z Na Z Nc e- EN ° ta e- EN ' ( - tl - ta) 



(T b ,p\N-0)\N a ,p) ^^ e -^* 2e -^(^ 1 )(-l)(* 1 +^cT C + 



Z N >Z N >e "o-'e 



Z Wb Z N ,e- E ^ 



E N I *2 e ~ E N' (tl-t2)( 



-E N ,(T-t 1 ) 



N a (p) 

N' a ,p M-p) K,6 
w b ,-p\NY(6) W a j 

W c ,6\N a (p) TT b ,- P 

N a ,p\n(-p) \N' C ,6) ^Z N ,Z Nc e- E K ta e- E ^- ta \-l) ta + 



-I) 



(-l)*^c + 



') ^Z JV ,e-^( tl - ta )e-^ (T - tl) (-l)( tl+ta )a c + 



( K,p\M-P) \N e ,0) ylz N ,Z Nc e- E ^e- E K^\-lf^). 

W = CttC -1 denotes the pion operator transformed by lattice charge conjugation and is 
the corresponding lattice charge conjugation number of the pion. N' is the parity partner 
state of the nucleon N. a c denotes the eigenvalue of the baryon state under a twofold 
application of the lattice charge conjugation. 

To interpret the matrix elements in the above formula we apply the LSZ formalism 
and evaluate the Green's function by employing the effective interaction Lagrangian of 
eq.(18), 



N a ,p 



Nc,0) = {9ic\ bc + g 2 cl bc )K(p)JZ Wh P((E Na - E Nc ),E Vl ), (22) 



and similar expressions for the other matrix elements. P(E,uj) denotes again the free 
particle lattice propagator eq.(7) and K(p) = |p|/(4rajv) is the appropriate kinematical 
factor. 

The function was fitted to the numerical data using a x 2 fitting procedure. Masses 
and energies were taken as input parameters from the 2-point correlation function analysis 
so that we have 4-fit parameters: two parameters proportional to the coupling constants 
9nnw an( i Qn'N'w an( i two proportional to the mixed matrix elements where the nucleon as 
well as the N' state contribute. An example for one of the fits is shown in fig. 2. 

The output of the fit is a product of Z amplitudes, the values of which are also taken 
from 2-point functions, CG coefficients and coupling constants. To calculate the CG 
coefficients c 1 abc and c^ Jc , we need to assign SU(4:)f quantum numbers to the states created 
by the baryon and the pion operators. For the 7r meson we take over the assignment of 
flavour quantum numbers discussed in the previous section. For the baryon, again using 
the staggered quark fields q, we can find an expression which enables us to investigate the 
quark content of the operator, 

N(y) = ^e ijk (?Aqi(yp A ) X {q 3 (y)[CT A ® C*T A ]q k (y)} . (23) 
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Here C denotes the lattice charge conjugation operator and a summation over colour indices 
k is understood. By transforming the operator N(y) to momentum space and applying 
the standard trace theorems for 7 matrices, it can easily be shown that this operator and 
the local baryon operator defined in [3] are the same. To find the correspondence between 
the lattice operator and members of the continuum 20-plet of baryonic | + states, the 
flavour content of the lattice baryon, which can be read off eq.(23), must be compared 
with the flavour content of the continuum baryons. This procedure is rather tedious for 
four quark flavours, leading to the result that the lattice baryon operator corresponds to 
a linear combination of three continuum states 

N Lau ^ J=(V3(A; + A° 2 ) - V2( Pl + p 2 ) - (St + St)). (24) 



The notations for the particle states are taken from ref.[8]. Since in SU(4:) F the 20- 
dimensional baryonic representation appears twice in the reduction of the representation 
of a 3-quark state we indicated by lower indices to which of these representations the 
baryon belongs. Specifically, S + lies in the 6-dimensional representation of the u,d,s 
SU(3) F subgroup of SU(A) F , in a 2-dimensional isospin SU(2) F representation and has 
isospin I 3 = +1/2. A belongs to the 3*-dimensional representation of SU(3) F , to a 2- 
dimensional isospin SU(2) F representation with isospin I 3 = —1/2. p denotes the proton. 
With the flavour content of the local pion already given in the previous section, the SU(4:) F 



CG coefficients c 1 abc and c 2 abc are found to be 



36^13 
1 1 



= (25) 



r 2 

^abc 



6V13 3V39 



To relate the coupling constant of interest §nn-k to the SU(A) F couplings §i and g 2 , 
we first consider the case (in the continuum) of a member of the baryon 20-plet belonging 
to a 2- and a pion state belonging to a 3-dimensional representation of the same SU(2) F 
subgroup. Then the strength of the baryon coupling is parameterised by gjwnC^ where 
c^ 2 ) denotes the corresponding CG coefficient in SU(2) F for this process. On the other 
hand, both the baryon and the meson are members of SU(4:) F multiplets and we can 
alternatively parameterise the coupling by g\C^ + according to eq.(18). Since these 

are just different prescriptions for one and the same physical coupling, we can write 

9nn-kC {2) = + # 2 C2 4) - (26) 

By fitting the numerical data in our simulation to the parameterisation function we 
extract a number g which is, according to eq.(18), given by g = g\C 1 abc + g-i^hc with known 
SU(4:) F CG coefficients. Using eq.(26) we can relate g to g^Nw through 

1 c{ 4) + c ( 2 4) g 2 /g 1 
9nn-k — S'tt^t— ; ; — 5 / — ' \ z ') 



c (2) c\ hc + cl bc g 2 /g 
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where everything is known except for the ratio g 2 jg\. To determine this number we used 
two different ways. 

First, we related the ratio to the experimentally known F/D ratio gpjgj^ in SU(3)f- 
This experimental value corresponds to the broken SU(3)f symmetry between u, d and s 
quarks in the continuum. However, this relation can be regarded as an excellent approx- 
imation for the symmetric case. We therefore conclude that the relation between gpjgj^ 
and g-ijgi can be based on the exact SU(4:)f symmetry. 

The relation between g^/g^ and g 2 jg\ can be found by a similar line of arguing which 
led already to the relation between g^Nw and 9i> 9'i- We choose baryonic and pionic states 
and regard them as members of 8-dimensional SU(3)f representations as well as members 
of an SU(4:)f 20-plet (baryon) and a 15-plet (pion). The baryon-pion coupling can now 
be formulated both in terms of SU(3)f and SU(4:)f couplings and CG coefficients. As 
the physical result is the same we have 

(4 3) ^/^ + cf)gi = (4 4 W<h + c[ 4) ) gi , (28) 

where c^ 3 ) denote the SU(3)f and the SU(A) F CG coefficients. However, we still need 
the ratio g\jg^ to find the desired relation. Therefore we consider a specific process, e.g. 
the coupling AA7T where one finds by explicit calculation that = as well as c!, 4 " 1 = 0. 
From this we can read off directly the relation 

„ _ „A "(3) / ~(4) _ A I ( 20 11 
9i-9d c d/ c i -9dI\ g g 

where the expression in brackets denotes an SU(3) singlet factor [8]. Inserting this result 
into eq.(28) we find g 2 /gi = 1.025 from the experimental value of 9fI9d = 0.8 (taken 
from [9], normalized according to de Swarts phase conventions for SU(3) Clebsch-Gordan 
coefficients [10]). This in turn enables us to find the desired relation between g and gNNw- 

g NNw = 1.6891st. (30) 

Based on these ingredients, we extract a value for g^Nw from our numerical data. We 

find 

g NNw = 13.8 ±5.8, (31) 

which is to be compared with the experimental value of 13.13±0.07. The result is consistent 
with the value quoted in ref.[5] gNNw = 12.7 ± 2.4 which was obtained by using Wilson 
fermions. 

The second method is based solely on lattice calculations. For this purpose we nu- 
merically analysed a second correlation function using the same baryonic operator but a 
different pion it'. This leads to a set of different SU(4:)f CG coefficients c 1 abc and c 2 abc and 
a different coupling g' = g\C 1 abc + g-i(? a hc so that we get two independent relations for the 
SU(4:)f couplings. Under the assumption of flavour symmetry gi and g 2 have the same 
value for all pion states in the 15-plet and we are able to solve for gi and g 2 separately. 




9d 



4^6' 



(29) 
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We used the spatially non-local pion operator [No. 14] for this purpose which corresponds 
to the continuum operator 



7T 




(32) 



from which we derive the CG coefficients 



abc 



61(2 + ^3) 
576^13 ' 
5(3 + 2^3) 



(33) 



abc 



18^13 



If we use the values for g and g' extracted from our numerical data, we are able to 
compute g^ldi, 



The large uncertainty of g-ij 9\ within the second method is mainly due to statistical noise 
in the numerical data. A possible explanation for the poorer quality of the NNir' 3- 
point function compared to the NNir correlation function, is the non-local character of 
the 7r' operator which requires the insertion of gauge fields in order to ensure local gauge 
invariance. 

The result for the ratio quoted above is smaller than the value g-ij 9\ = 1-025 we 
obtained previously, though both are compatible within the errors. But we have also to 
bear in mind that the second value was obtained under the assumption of SU(4:)f flavour 
symmetry on the lattice. However, the discrepancy between the masses of the lattice 
Goldstone pion and e.g. the mass of it' we found in our spectrum analysis indicates that 
SU(4:)f flavour symmetry is not yet completely restored for our lattice parameters. In 
order to roughly estimate the effect of flavour symmetry breaking, one could argue that, as 
a result of the higher mass of the it' state, the virtual transition N — ► N + ir' is suppressed, 
thus leading to a smaller value for g'. If g' increases, the ratio g-ij 9\ rises and comes closer 
to 1.025. 
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Figure 1: The numerical data for the 3-point function in eq.(2) with the local p operator 
(positive values are marked with crosses, negative values with squares). A fit from t x = 1 
to ti = 23 without ti = 4 (the time slice of the pion source) is shown by the full curve. 
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Figure 2: The numerical data for the 3-point function defined in eq.(19) with the lo- 
cal Goldstone 7r operator (positive values are marked with crosses, negative values with 
squares). Also shown is a fit from t 2 = 1 to t 2 = 23 without t 2 = 8 (the time slice of the 
nucleon). 
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